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Abstract 

We discuss a systematic way to dimensionally regularize divergent sums arising in field 
theories with an arbitrary number of physical compact dimensions or finite temperature. 
The method preserves the same symmetries of the action as the conventional dimensional 
regularization and allows an easy separation of the regulated divergence from the finite 
term that depends on the compactification radius (temperature). 


1 The problem 

In a variety of problems one has to deal with formally divergent sums, usually related to 
Feynman diagrams with one or more discrete momenta, as in the case of theories with 
compact extra-dimension, finite-size scaling theory in critical phenomena, thermal field 
theory. It is crucial to find a regulator that preserves the symmetries of the problem 
and leads to a simple computational procedure. When the momenta are not discrete but 
continuous, such a procedure exists and it is the well-known dimensional regularization 
(DR) of integrals In this paper we discuss a systematic way to obtain the dimensional 
regularization of an important class of sums, following closely the analogy with the case of 
a continuous variable |p. Even if practical recipes to deal with particular examples have 
been given in the past and the use of special functions to this purpose is not new 
our aim is to provide a general method, which extends and in a sense justifies the analysis 
of ref. 1^, where the idea of dimensionally regularized series was applied to a number of 
different loop sums. In |]l^] the approach to dimensional continuation was adopted and 
combined with complex analysis techniques to delineate a general procedure, restricting 
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however to the case of only one physical compact dimension. In |]^ the use of the 
Mellin transform, together with dimensional regularization was introduced to deal with 
asymptotic expansion of series in thermal field theory, formally the same problem as 
discussed in |n|. 

We propose to define sums in complex dimension using the analytic properties of a 
generalized zeta function, resulting in a simple method where the regulated divergence 
can be easily separated from the finite part. Our technique has some common points with 
the well known zeta-function regularization |^, which leads to quite similar calculations 
but is well distant from the spirit of analytical continuation in the number of dimensions. 

We hope that this work may contribute to clarify some debated aspects of regular¬ 
ization in extra-dimensional models, where it is crucial to preserve the symmetries of the 
action |13-p^- In |]^| it was suggested that the finite result obtained by Barbieri, Hall 
and Nomura |]^| for the radiative correction to the Higgs mass coming from the Yukawa 
sector of their model was a regularization artifact. Unfortunately, the authors of ref. |14] 
made use of a sharp cut-off on the series, which explicitly breaks the supersymmetry of 
the model of ref. []^| and invalidates their argument. The result of |]^] was reobtained 
in 1151 by using a thick brane as a regulator for the series and in ||^ in two different ways, 
with a Pauli-Villars and by using dimensional regularization. In this work we give a 
formal procedure to introduce the dimensional regularization of a series and demonstrate 
in detail some important properties. In particular, we show that there is no ambiguity in 
exchanging the series and the integral over loop momenta if both are properly regularized. 

The paper is organized as follows. In Section we briefly review standard dimensional 
regularization of integrals and give our rules for extending it to the case of a series by using 
a generalized zeta function. We discuss the infinite radius limit in Section showing that 
the (regulated) divergence equals that of the corresponding integral and as such does not 
depend on the radius. A representation using theta-functions is also given in Section |3.2| , 
which is useful to perform explicit computations of sums in a class of physical problems. In 
Section ^ we show how the method works on an explicit example, computing the Casimir 
energy of a massive scalar field. In Section we draw our conclusions. 

Finally, in the appendices we briefly recall the definition of the Mellin transform ( 0 ), 
present the analytic continuation of the generalized zeta function, together with its asymp¬ 
totic behaviour and discuss some important properties of the dimensionally regularized 
series 


2 Dimensionally regularized series 

Dimensional regularization of integrals was introduced in Q as a simple tool to manage 
the divergences that arise in (perturbative) field theory, preserving the gauge symmetry. It 
was later derived as an axiomatic procedure by Wilson |||], see for a detailed discussion. 
The problem is to give a meaning to the integration of a function f{p^) over a space of 
complex dimension d, getting the usual result whenever d is an integer and the integral 
exists in the ordinary sense. Without going into too much detail, one can define the 
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integral through the formula 


/ TOO 27 r“/^ 

d'^pfiP^) = dpp'^-^ f{p^) ; VLd = Y[d/2) ’ 

for all (complex) values of d for which the integral converges and then analytically con¬ 
tinue the result to the desired value. Typically one has to cope with ultraviolet (UV) 
divergences, which are cured by considering a sufficiently small (real part of) d and ap¬ 
pear as poles in the final expression for the integral. It may happen that the value of d 
that makes (|^) UV-convergent is so small that the integral diverges in the infrared (IR): 
in this case one dehnes the integral subtracting the leading behaviour of the function / 
for p^ —> 0. Let us suppose that f{t) has the following asymptotic expansion for t —> 0: 

/(t) = aot“° + ait“i+02t“^ + --- 

(^) 

ctQ < cti < q;2 < • • • 


and f{t) ~ t ^ for t —> oo. Then, for —ctg < Re d/2 < min(—a^-i, p), q a positive integer. 


the integral is defined by 



/* roo 

jdSf{p^)=^d dpp^^-^ 

1 - 1 

to 

9-1 

-^ak{p^T^ 

k=0 


If —ao < p subtractions are not really needed and the integral on the r.h.s. of eq. ® 
converges for — ao < Re d/2 < p; in this case eq. (|3[) is simply the correct analytic 
extension to the interval —Og < Re d/2 < —ag-i, q > 0. On the contrary, when —ao > p, 
there is no value of d that makes (|^) both UV- and IR-convergent and (Q) becomes a 
dehnition. 

This general procedure does not work in the case of sums, given that there is no closed 
form for what corresponds to the “solid angle” 11^ in the case of a hypercubic lattice and 
one cannot directly take advantage of the spherical symmetry of the function /. We then 
proceed as described below, following closely the analysis of |^] for the integrals. 

We want to give a meaning to the expression 

neJA 


where the prime means that n = 0 is omitted in the sum and / is assumed to be continuous. 
We restrict for simplicity to the case of a scalar function (the case of tensorial functions 
can be addressed in the same way as for standard DR of integrals) and require covariance 
of our result under (discrete) rotations of the hypercubic lattice. In the case of only one 
variable, this implies that the function / depends only on the norm of the vector n and 
not on its direction. The general case in which / depends also on external momenta is 
considered below. 
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An operation of summation in arbitrary complex dimension is uniquely determined [| 
by requiring the following basic properties, valid for standard summation: 

1. Linearity: for any complex numbers a,b 

'^'[a f{n)+b g{n)]=a'^'f{n)+b'^'g{n) ; (5) 


2. Invariance under lattiee translations: for any vector q 

y^/(n + g) = y]/(n) . (6) 

riEZ^ 


Notice that the scaling axiom ^ required for dimensionally regularized integrals does 
not hold for the series. As for the standard case of integrals, vectors are thought to he in 
an infinite dimensional space, with the difference that now each component of the vector 
has an integer value. The dimensionality d is introduced by the sum operation with the 
requirement that if d is a positive integer all vectors collapse in a d-dimensional subspace. 
When the function / depends also on some external momenta g* (actually only through 
the scalar products (n • qi), {qi ■ qj) being a scalar function), one can proceed again in 
complete analogy with Q: it is always possible to find an A^-dimensional sublattice 
(with N finite being the external vectors in finite number), which contains all the external 
vectors. Let us decompose n into a longitudinal and a transverse part with respect to Z-^: 
n = n|| + n^, so that 

'I2f{n,qi)= y] y] f{n'i + nf^,n\\qi,qi-qj) . (7) 

nEZ^ riji EZ-^ nj_ EZ^“-^ 

The outer sum in the r.h.s. of (|^ is a standard series on the lattice Z^, while the inner 
one can be defined through (^), / now being independent of the direction of n±. Then 
there is no loss of generality in reducing to the case in which / is a function of only the 
dummy variable. 

Let us now describe an explicit procedure for summing in complex dimensions. Given 
a function f{t) continuous for t G M+, let us assume that 9{t — a)f{t) (where 9 is the step 

^Except for an arbitrary normalization, which can be fixed on a set of basis functions |^,^. We 
require the usual result of integer dimensions 

= d^(s) , 


to hold for all d G C (see the appendix for the definition of the d-function). 
^It states that V s > 0, 
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function) is Mellin-transformable with fixed 0 < o < 1 (see Appendix^ and Ref. |p^ , p^ 
for details on the Mellin transform), i.e. 3 p G M such that 

< oo ya < p . 

Then the Mellin transform 

poo 

Ma[f,s]=M[e{t-a)f{t),s]= (8) 

J a 

is an analytic function for s G C if Re s < p. In particular, the assumption of the 
existence of the Mellin transform excludes from our discussion those functions f{t) growing 
exponentially for t +oo. The inversion theorem guarantees that 

-j pc+ioo 

e{t-a)f{t) = — / ds Ma[f,s] . (9) 

^ Jc—ioo,c<p 

Usually, the Mellin transform of the function / is defined in the strip —ao < Re s < p, if 
fit) fOLQ t —> 0; choosing the parameter a > 0 allows the lower limit —ao to be sent 
to —oo. Let us now assume that there exists an integer value of d such that (§) exists in 
the ordinary sense. It is easy to realize that if it is the case, then d/2 < p and 





1 

2tt i 



1 

(^2)s 


( 10 ) 


where exchanging the series with the integral is allowed in the domain of uniform conver¬ 
gence (see below). We define a generalized ((-function (a particular Epstein’s zeta function, 
see and |p^ ) 

• ( 11 ) 

whose properties are discussed in Appendix p. This series converges uniformly in any 
closed subset of the line Re s > d/2, i.e. the function Cis,d) is analytic in the half-plane 
Re s > d/2. Thus we may give the following representation of our original sum; 

^ [ds C{s,d)Ma[f,s] ; T = {Re s = c-,d/2 < c < p} . (12) 


Let us observe that both M.a[f,s] and Cis,d) can be analytically continued outside their 
definition domains. Their analytic continuations will generally have singularities in the 
complex plane. Now the recipe to define the sum (Q) even for d/2 > p is clear: it is 
simply the continuation of integral ( 0 . where we consider the analytic continuation of 
Aia[f,s] and Cis,d) in the integrand. In the same way we can define (^) by using the 
representation (|^ even if there is no integer value of d such that d/2 < p (for example 
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if f{x) ~ x“ for X ^ oo and a > — 1) and the same considerations apply for complex d if 
one considers Re d instead of d in previous relations. 


When Re d is increased towards values greater than p, the integral (^) gets a residue 
contribution from the pole of Cis,d) in s = d/2: as we will show below, this term is 
divergent for Re d/2 > p and coincides with the infinite radius limit The remaining 
complex integral along the contour T = {Re s = c;c < p < d/2} corresponds to the 
finite radius-dependent part; in the following section we give a practical recipe to evaluate 
it using the r^-function. Note that the operation (12) respects the required properties: 
linearity follows from linearity of the Mellin transform; translational invariance follows 
from translational invariance of the usual sum over the transverse subspace, if this is 
taken sufficiently large to contain the vector q in eq. (P). In Appendix ^ we prove some 
remarkable properties which hold for the sum over complex dimensions. 

Finally, a comment on ([T^) is in order to clarify the meaning of the parameter a: 
as already said after eq. ( 0 ), going to too small d without performing an appropriate 
subtraction in the case of the integral would introduce a spurious IR divergence which 
has no physical meaning. In the same way, had we not introduced the parameter a > 0, 
the Mellin transform 


fit) 

would have been divergent (for any s) in all cases in which the integral corresponding to 
the series needs IR subtraction to be dehned in dimensional regularization. Setting o > 0 
is a simple way of avoiding IR subtractions for the series. We prove below that the final 
result does not depend on a, as one expects having the initial series (Q) no IR problems 
at all. 



3 An alternative representation 


Once given the definition of a series in complex dimensions, we present an alternative 
representation of (12) to show that 

- the result is independent of the actual value of 0 < a < 1; 


" the UV dimensional poles are the same as those of the corresponding integral; 


- the property (analogous to that for ordinary sums) 


lim 

R^oo 





is formally valid. 


(13) 


The final formula will be easier to handle for a numerical evaluation or an expansion in 
the parameters. 

^If the physical value of d is less than 2p the series is convergent and there is no need to 
(dimensionally) regularize it; in this case the representation ( p^ is still useful for an explicit 
evaluation. 
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Figure 1: Position of the contour P and Ph the poles for Re s > p, are those of the 
continuation of M.a[f, s]. 


3.1 The infinite radius limit 


Our strategy is to obtain first a compact expression for Re d/2 < p and then analyti¬ 
cally continue this result to the case Re d/2 > p. Introducing the radius R explicitly, 
definition (|l^) becomes (0 < a < 1): 


Sid,R) = ^ = ^J ds as,d) R^^-^Ma/R2[f,s] 


(14) 


where the contour P has to be fixed according to (p^. Let us assume that f{t) has the 
asymptotic expansion (§) for t ^ 0 and f{t) ~ t~^ for t —> oo, even if the derivation goes 
the same way for more general expansions. We can take 


P = {Re s = c] —Oq < Re d/2 < c < min(—ag_i,p)} , 

q being the smallest integer such that —Og < Re d/2 < p (see Fig. |T|). Dehning^ 

9-1 

/sub(i) = f{t) - '^ak , (15) 

fc =0 


it follows, for —Uq < Re s < min(—ag_i,p): 

9—1 / n \ 

Ma/R2[f,s]=M[fsnh,s]- df/sub(i)t^“ “ X] “T- [^) ' 

Jo s + ak 

Note that the Mellin transform in the first term is IR-convergent thanks to the subtraction 
made on / and the second term is analytic for Re s > —aq. Using the asymptotic limit 

■^For q = 0 the subtraction terms are omitted and we mean min(—Q!_i,p) = p. 
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C{s,d) —> 2d for Re s ^ +oo (see Appendix |^, and moving the contour T at infinity in 
the positive half-plane, one easily finds that the integral in the complex variable s of the 
second term vanishes, since a < 1. Integrating the third term along F, again closing the 
contour in the positive half-plane, one gets instead a residue contribution from the poles 
s = —Ofc, k = 0,..., {q — 1). We end up with 




(17) 


This expression represents a definition of the series for —Og < Re d/2 < min(—cog-i, p) 
alternative to (0). where the use of the IR cut-off a has been replaced by an appropriate 
subtraction of the initial function, much in the same way as one does in standard DR. 

If — ao < P the IR subtractions are not really required to define the series: choosing 
d in the strip —ao < Re d/2 < p {q = 0), eq. (pA|) coincides with (1^) with a = 0. In 
this case eq. (0) with q > 0 represents the correct analytical extension of the series to 
the interval —a^ < Re d/2 < —a^-i, but it becomes non-trivial when —ao > p and the 
subtractions are really necessary to give a meaning to the expression. It is interesting 
to notice that Ad [/sub, s] in the first term is nothing but the dimensionally regularized 
integral of / in d = 2s dimensions, up to a solid-angle factor. 

The limit R ^ oo cannot be extracted immediately from (^), as one must first 
extrapolate the result to values of Re d/2 greater than —ao- To analytically continue (ITtI) 
to Re d/2 > min(—ag_i,p) we have to take into account that increasing Re d/2, the pole 
s = d/2 of C('S, d) crosses the contour T and gives a residue term. To simplify the analytic 
continuation, we can rewrite (|^ in a suitable form, evaluating this residue term from 
the very beginning, i.e. moving the contour F into a new one (see Fig. [T|) 


F' = {Re s = c; —Og < c < Re d/2 < min(—a^-i,p)} 


so that 


jr-d./2 1 f 

Sid,R) =f(^-^[/sub,d/2] + — ds as,d) + 


q-l 

+ E 

k=0 




J^2o'k-\-d 


C(-afc,d) 


(18) 


The continuation of this expression can be easily found by considering those of M. [/sub, d/2] 
in the first term and of C('S,d) in the second and third ones, while the contour is always 
given by F' = (Re s = c; —ag < c < min(—ag_i, p. Re d/2)}. 

Dimensional poles of this expression may arise from the first and third terms on 
r.h.s.: the former is i?-independent and exactly corresponds to the integral (^ defined 
in dimensional regularization with the proper subtraction. Its possible IR divergences 
are cancelled by the poles of the last term for d = —2ak, k = 0, - ■ ■ , (g — 1), which may 
be “physically” accessible, i.e. be positive, if some a^ is negative. The final result is 


thus IR-finite, as it should be, the series being free from IR divergences. From (18) it is 








evident that the UV divergences are the same for the series and the corresponding integral. 


Finally, when Re d > — 2ao the infinite radius limit can be safely extracted from (^) 


leading to eq. (13) as promised. 


3.2 Evaluating the finite part with a i9-function 


A general formula for the complex integral in ( [iq ) can be obtained under particular 
assumptions. Performing the appropriate IR subtraction, one always reduces oneself to 
the computation of the integral 


1 

27rz 


dsC{s,d)R^^-’^M[f,s] , 


(19) 


where, if the function / has the usual asymptotic expansion (|^ and ^ for t —> oo, 

then r = {Re s = c; —ao < c < min(/?,Re d/2)}. Let us suppose that 


1. f{s) is analytical in the half-plane Re s > c, c < 0, 

i.e. /O > 0, 

pc-\-ioo 


2. lim f{s) = 0 

Re s—>+oo 


3. 


1 

27r i Jc—i 


ds f{s) < oo for t > 0. 


c —200, C>C 


Under these hypotheses the inversion theorem (see |jl^) guarantees that the Laplace anti¬ 
transform 


1 1 * 0+200 
= — / ds f is) 

ZTT I Jc—ioo.c>c 


( 20 ) 


exists for t > 0 and 


POO 

f{s)= / dt C~^[f,t] Re s > c . (21) 

Jo 

Using these results, the Mellin transform can be expressed as follows: 

POO POO POO POO 

■M[f,s]= dtR~^f{t) = j dtf~^ I dye~^yC-^[f,y]=T{s) I dy £"^[/, y] y"* . 

Jo Jo Jo Jo 

( 22 ) 

This relation holds when —oq < Re s < p and it is therefore valid along the contour 
r of eq. (^9|). If ao > 0, Re s can be negative and = 0, for A: = 0, • • • , [aoj, 

[xj = maxjn G N | n < x}. Then the simple poles in s = 0, —1, —2, • • • , — [aoJ of r(s) 
in (|^) are cancelled by corresponding first order zeros of the integral of the Laplace 
anti-transform 


dyC ^[f,y]y^ = {-!)'" = 0 for A: = 0, l, • • • , [aoj if ao > 0 
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resulting in a Mellin transform well defined in the whole interval —ao < Re s < p. 
Using and the analytic extension of C(s, d) in terms of the function (see Appendix]^) 


Cis,d) = 


TT” 


r(s) ( s — d /2 




d\y) - 1 


it is not difficult to compute the complex integral in (|I^) by applying standard residue 
techniques. The result is: 


1 

27ri 


T [ dsCis,d)R^-^-'^M[f,s] = 

I Jr 

/ OC 

dy [/, yTiR^] [/, TiR^/y] y~ 

/ oo , 

dy ['d'^(y) - 1] {c~\f,TTR^y] +y'^/‘^~‘^C-\f,TTR 


TT 

IT 

-^2 


(23) 


The two integrals converge under the assumed hypotheses Equation (|^) is a useful 
representation for the complex integral in (|l8|) , which allows an expansion of the result in 
terms of the parameters in / or even to perform a numerical integration. The results of 

It is worth noting that one 


ref. I IT I are easily obtained as particular applications of 
could have derived the same result by giving a definition of the series directly in terms 
of the function, never using ((^-functions, and this is exactly the procedure adopted in 
ref. II^. Of course, to define the series using functions one has to assume the same 
hypotheses we imposed on the function / to guarantee the existence of the Laplace anti¬ 
transform. In this respect, our definition seems more general and more suitable to isolate 
the divergent from the finite term. 


4 A sample computation: the Casimir energy 

To give an explicit example of our procedure, we compute the Casimir energy £{R) for 
a scalar field with mass m and periodicity conditions in a space geometry T'^ x M^, 
corresponding to d compactified dimensions with radius R. This is defined as the R- 
dependent part of the zero point energy V{R) of the (Z)-|-(i)-dimensional theory multiplied 
by the volume {2'kRY of the compact space: 

£{R) = {27TRf [V{R) - U(oo)] , (24) 

®This is easily understood from the asymptotic behaviour of the Laplace anti-transform 
(see 0): 

{ fP ^ f —> f) f f P f —^ r)Q 

, . , which follows if /(t) ~ < 

t-(“0 + l) t^oo t^O 
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Continuing the physical dimensions D, d to generic values 

D = D — e] d = d — rj , 


( 25 ) 


the one-loop contribution to the zero point energy is given by 

= 5 (Sflp E / ‘oe {/ + ..VR" + m") ^ (26) 

ndlA 

To extract £{R) we first compute V{R), then extract the limit e —> 0, r/ —> 0 and finally 
multiply by the volume factor. 

For convenience let us call Vo{R) the zero mode contribution in (^). We will show 
explicitly that performing first the series and then the integral or vice versa leads to the 
same result (see property 2 in Appendix and there is no ambiguity. Let us perform the 
integral first. The function f{x) = log(x + c) has the asymptotic behaviour f{x) ~ logx 
for X —> oo and f{x) ~ const, for x ^ 0; this means that the dimensionally regulated 
integral needs subtractions to be defined. Performing one subtraction, i.e. considering 
/sub(^) = log(x/c + 1), we can fix —2 < D < —1, obtaining 

V{R) - Vo{R) = -- f'+’i ^r(-D/2) ■ (27) 

ndlA 


Because D is negative and the function g{x) = (x + m^)^/^ has the asymptotic behaviour 
g{x) ~ x-^/^ for X —> oo and g{x) ~ const, for x ^ 0, the series is dehned by (^) with 
r = {Res = c;0 < d/2 < c < —D/2 < 1} and without any subtraction: 


1 


V{R)-Vo{R) = --ix 


TT 


2 ^ (27r 


DI2 

15n-D/2) 


1 


1 


ds C(s, d) R?^M[g, s 


{2ttR)‘^ 2m 

Even if the Mellin transform is easy to derive 

r(s)r(-i:»/2 - s) ^ 2 s . d / 2 +s 


(28) 


M[g,s\ = 


n-D/2) 


( 2\^ 

[m j 


the complex integral along the contour T cannot be solved with simple residue techniques 
because of the non-trivial behaviour of the integrand at infinity. However, g(s) is analytic 
in the half-plane Re s > —m? and the requirements needed to apply the representation in 
terms of the -d-function are fulfilled. From eqs. (p^) and (^) with a Laplace anti-transform 


C- \g,y] = 


Q-in.^y 

^{-D/2) 


-D/2-1 
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we obtain 


V{R) - Vo{R) 


1 , / 

-_r I 

2^ {2TT)D+d^ 


D + d\ 

2 y 


(^2){i^+rf)/2 


+ 


1 /i"+^ 

2 (27ri?)^+'^ 


/ OO 

dy [^'"(y) - 1] 


y-(D+d)/2-l 
f^--Ky{mRf y 


-D/2-1 _|_ lyy{D+d)l2-l 



We recognize the infinite radius contribution in the first i?-independent term, while the 
remaining terms are finite. More exactly, the integral in the second line is convergent 
when Re D < 0 as supposed from the beginning to regularize the dimensional integral. 
Therefore, when Re D is increased to (physical) values Re D > 0 a divergence appears 
which, however, depends only on D but not on the value of d. This means that in some 
way this must be a “zero mode” divergence and this becomes evident by rewriting 


^ = [7r(mR)2]^'^^r(-T)/2) - ^ 


- / . (30) 


The first term on the r.h.s. of (|3b| ) exactly cancels the zero mode contribution Vo{R) 
in (29) and we get the final result 

V(R) = f-:^) (m2)(o+<i)/21 X 


(27r)'^+'^ 


2 (27rR)^+'^ 


dy fy-{D+d)/2-i ^y-D/2-i\ _ 


i 

dy [^'^{y) — 1] y~^/VyiD+d)/'^-^'^ | 


(31) 


The first term is the ordinary divergent renormalization of the cosmological constant, 
which can be put to zero with a suitable counterterm if we accept the usual fine tuning. 
Whatever scheme of renormalization one chooses, the radius dependent part V(R) — V(oo) 
of the zero point energy is non ambiguous and finite. We can therefore safely extract the 
limit e —> 0, 7/ —> 0 and insert the result for V(R) — R(oo) in eq. d^) to extract the 
Casimir energy 


T(R) = 


2 (27rR) 


D 


dy fy-{D+d)/2-i ^y-D/2-i\ _ 


/ CXD _ . 

dy [d'^iy) - 1 ] + f,-^{mRf/yy{D+d)/ 2 -i 


(32) 


The same result can be obtained by performing first the series. Again, subtractions 
are necessary if one does not introduce an IR cut-off in the Mellin transform. The function 
/(x) = log(x + + m^) has an expansion around x = 0 as in (|^ with oq = 0, ai = 1; 
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we can therefore apply definition ( 0 ) with one subtraction and, choosing the contour 
r = {Re s = c; —1 < (i/2 < c < —1/2}, we obtain: 


V(R) - V„(R) -- 


+ 


C(0, d) log + w?) + 

2(2.Kf J {2.)0 2h X* 


1 


/ 


(27r)^ 
d^p 


(33) 


with 


Ai[/sub,s] = - [p^+ m^]^r(s)r(-s) . 

The complex integral along the contour T cannot be solved with simple residue techniques 
nor we can apply the li-function representation (/sub diverges at infinity), but we can make 
use of a standard trick. Namely, defining 

S{p^ + ^ C(s, d) R^*Al[/sub, s] 

we know how to sum the series dS{p^ + m?‘)/dp^ because dfsnh/dp^ = l/(x + p^ + m?) 
goes to zero for x —> oo. Then we can deduce the expression of S{p^ + w?) except for 
an unknown function independent of p^, which however is irrelevant as its dimensional 
integral in p gives zero. In particular we have: 


1 

27ri 


d 


ds C,{s,d) R^^^Al[/sub,s] 


d 


1 


(34) 


(is C('S,d) R -r^A4[/sub,s] - C(0,d)-2 , 2 

dp^ pz _|_ 


where the contour T has been moved into a new one, T' = {Re s = c; d/2 < 0 < c < 1} 
and a residue contribution in s = 0 has been isolated. The first term corresponds to 
the (dimensionally regularized) series of the function dlog{p^ + n?/R^ + rn?‘)/dp^, defined 
without subtractions by choosing the contour T', and it can be easily computed using the 
function representation (p^ . We do not show the details of this calculation but we only 


notice that using (^) to extract S{p^ + m^) and finally plugging the result into (^), the 
residue term in the former equation cancels the subtraction term of the latter and we end 
up with 


V{R) - Vo{R) 


d^p 

{2tt)^ 


IT 


DI2 


T{-D/2) [p2 + m2]'^/^ + 


(27r) 


{2ttRY 


/OO y . 

/•OO . 

I dy [d'^iy) - 1 ] + ^-^{m^+p^)RVyydl2-1 


(35) 


The result of the dimensional integration over p coincides with eq. (29), from which the 
final expression (31) follows. 
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5 Conclusions 


We discussed in detail a general procedure to dimensionally regularize divergent series. 
The novelty of the method consists in using a suitable combination of two well-known 
tools, i.e. Mellin transform and analytic extension of special functions, to provide a 
continuation of the series in the number of dimensions. 

The virtue of conventional dimensional regularization is to preserve all the symmetries 
of the action that do not depend on the dimensionality, in particular gauge invariance and 
supersymmetry (if the dimensional reduction scheme is used). The same happens with 
dimensional regularization of sums, making this technique a natural choice to handle di¬ 
vergences of field theory. In this respect, our analysis should contribute to clarify some 
controversial aspects about the computation of quantum corrections in supersymmetric 
theories with extra dimensions EHl- In particular, we have shown that there is no 
ambiguity in exchanging the series and the integral over loop momenta if both are con¬ 
sistently regularized with dimensional regularization (see Property 2 in Appendix |y and 
the example in Sec. The same is not true for instance if the sum over Kaluza-Klein 
modes is truncated by a raw cut-off. 

Our definition of sums in complex dimensions by using a generalized zeta function is 
particularly suited to isolate the divergence from the finite part and applies to a large class 
of functions. Moreover, it leads to simple computations and it is valid for an arbitrary 
number of physical compact dimensions. The idea of dimensionally regularized series was 
applied in ref. to compute a number of loop sums by using a representation in terms of 
??-functions. We obtain the results of ref. [l^ as particular cases of our general formulae. 
A different method was proposed in [^, which however applies only to the case of one 
physical compact dimension. Our procedure has no such limitation and it is therefore 
more general; in the case of only one physical compact dimension it gives the same result 


as |H], as we have explicitly checked for particular examples. 

Although the class of functions we considered is not the most general one, the basic 
idea may be applied to more complicated cases by using suitable zeta functions. In par¬ 
ticular, our formulae are specific to series that appear in theories with toroidal compact 
dimensions and scalar fields with simple periodical conditions. In the case of twisted 
periodic conditions or theories with fermions, our method should be easily extended in¬ 
troducing the following zeta function 


C(s,ci| a) = 


1 


[(n -I- a 


|21s 


-1/2 < a < 1/2 


whose properties are sketched in Appendix p. The way to treat spinors in the dimensional 
continuation is similar as in conventional DR, see 10 - 

Even when the compact manifold is not toroidal (see for example |^) or the single 
dimensions have different radii, we see no obstacle in principle to apply the approach of 
dimensional regularization, maybe taking advantage of more general special functions such 
as those introduced in |^. In these cases, however, computations may become involved 
and particular recipes specific to the case may be simpler to use. 
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A Mellin transform 


Let us briefly recall the definition of the Mellin transform MM- Given a function /, if 
a, /9 G M, a < /? exist such that 

\f{t)\ t^~^ < 00 ; \/p a < p < (5 , 

then one can define the Mellin transform of /(t), 

poo 

M[f,s\= I dtf{t)f~^ , (A.l) 

Jo 

which is an analytic function of s G C in the strip a < Re s < (5. This inversion formula 
holds 

1 pc-\-ioo 

f{t) = — ; a<c<(3 . (A.2) 

iTT X Jc—ioc 



Obviously it is a linear integral transformation, 
markable properties, see, for instance, |n|. 


It has very useful applications and re- 


B Properties of (^(s, fi) 


We define, for integer d 


C{s,d) = 


1 


(B.l) 


which is absolutely convergent for Re s > d/2. This function is a particular case of 
the more general zeta functions (see, for example, and [^). Let us write an explicit 
expression for the analytic continuation of C{s,d) to the whole complex plane in s and, 
eventually, extend its definition also for all complex d. The function (related to Jacobi’s 
^3 @), 

+ 00 

E ’ (b.2) 

n =—00 


15 





has the following modular property, which is easily derived from Poisson’s resummation 
formula [jT^: 




By means of 


we may write 


1 1 


{n^r r(s) J„ 


f 


dt ^ e 


TT 


s roo 


pc 

C(s, d) = / dt - 1 

r(s) Jo 


(B.3) 

(B.4) 

(B.5) 


incidentally, this tells us that the Mellin transform of — 1 ] is vr ® r(s) ("(s, d). 


For Re s > 0, Re s > d/2, it follows from (B.2) that 

r -1 r 

Jo ^ ^ Ji L 


+ 


1 1 


which gives 


C{s,d) = 


TT 


r(s) i s — d /2 


1 r 

-s^L 


dt (t^~^ _^^d/2-s-l 


s — d/2 s ’ 


- 1 


(B.6) 


This expression represents the analytic continuation of (/{s, d) in both s and d. It is easy 
to demonstrate the following properties; 

1. Ci^jd) is a meromorphic function, having a simple pole for s = d/2, d ^ 0, with 
residue 


Res {C('S, d); s = d/2} = 


TT 


d/2 


r(d/2) 


2. C(s,0) = 0, Vs E C; 

3. C(O)d) = —1, for d / 0; 

4. C{—n, —2n) = (—I)” nl/vr"', for n £ N, n / 0; 

5. C(“^) d) = 0, for d / —2n and n E N, n / 0; 

6. Given the symmetry s^d/2 — sof the terms inside the brackets in (|B.6| ), we have 

r(s)C(s,d) _ r(d/2-s)C(d/2-s,d) 

Vr^ yj-d/2—s 

This is the so-called reflection formula, which is also valid in the more general case 
of Epstein zeta functions Q . 
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We now derive a “convolution” property of used in Appendix ^ eq. (|C.3| ). For 

p,q > 0 and Re s > (p + q)/2, 


C{s,P + q) 


E'E' 

m£ZP n&i 

E'E 

meZP neZ'J 


1 

(m^ + 

1 

(m^ + 


S 


S 


+ E' 

meZP 


1 

{m?Y 


+ E' 

n£Zl 


1 

{v?y 


+ C(s,p) + C('S,g) 


(B.7) 


We can write the argument of the double sum using (B.4) and then introduce the 
Mellin representation for one of the two exponential factors, getting the result 


_ _ 1 1 noo _ _ -j pc-\-lOO 

E' E' („., „2). = fM / * *'■' E'«“ rw 

V ' / I J -/O m£'EP Jc—ioo,c>0 


2 j.\ —w 


mGljP nGlj^ 


1 1 

r(s) 27ri 


^c+zoo 

7 / (ilf; r(i(;) r(5 — tf;) C(s — g) 

^ Jc—ioo 


(B.8) 


where, in the last line, p/2 < c < Re s — q/2. Taking into account relation (BT), we 
have PI 


2 rc+ioQ 

— / dwT{w)C,{w,p)T{s-w)C,{s-w,q) 

ZTT % Jc—ioo 


(B.9) 


T{s)Cis,p + q) -T{s)Cis,p)-T{s)Cis,q) ■ 


This relation is still valid if one considers the analytic extensions of the functions involved. 

Next, we determine the asymptotic behaviour of C{s,d) for large Re s and d > 0. Let 
us observe that we may write 

OO 

= l + 2(ie-^‘ + ^ A4(d)e-“''^‘ , (B.IO) 

k=0 

where A4(d) are real coefficients and 4 < ao < ai < • • •. For large k, the series is 
asymptotic to the one with ak = k^ and Mk{d) is the number of points in lY with a 
distance from the origin bounded between k and {k + 1). So Mk{d) ~ k‘^~^ (for k ^ 1). 

For Re s > d/2 we may insert ( |B.1(]| ) into ( [B.5D , obtaining 

OO 

C{s,d) = 2d + Y,^k{d)a^^ , 

k=0 

SO that 

OO 

|C(s,d)-2d| ^ A4(d)(afc/ao)"^"" • (B.ll) 

k=0 

®Note that it is not possible to close the contour of integration at infinity, given the asymptotic 
behaviour of the integrand. 
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The series on the r.h.s. is always convergent for Re s > d/2, and therefore Cis,d) — 
2d ~ for Re s —> +oo. Using the reflection formula we immediately obtain the 

asymptotic behaviour for Re s —> —oo: in this case the function is unbounded [|: 

c(-..d) . (-i)W 

A useful generalization 

In some cases of interest (such as in theories with fermions) it is useful to introduce the 
following generalization of our zeta function. Given a constant —1/2 < a < 1/2, we define 


C(s,d| a) = 




[{n + a^Y 


(B.13) 


where now Y/' means that for a = 0 we omit the term n = 0 in the sum. This series 
converges for Re s > d/2. Obviously ({s, d | 0) = C(s, d). To get the analytic continuation 
of this function one proceeds as in the case of Cis,d). Introducing the d-function ^ 


+ 00 


'&{t I a) = 


(n+a)^ 


(B.15) 


from the Poisson resummation formula easily follows the modular property 


^—TT t 


The analogues of eqs. 


= —rjY 
and (|B.6|) are 


-'&{l/t\iat) , 


C(s,d| a) = 


TT 


r(s; 


f 


dtt 


5—1 


I a) - (5a ,0 


C(s,d| a) = 


TT 


r(s) \ s — d/2 s 


-(5a ,0 + / dtt 


i: 


.5—1 


rit I a) - (5a ,0 


+ 


^e-^<^^d/t^<i(t\ia/t) -ij I 


’’We use the fact that for x S R+, m= \x\, q = x — [xj, one has 

^T{q+l)n-q) 


n-x) = {-ly 


r(x + i) 


T[-x) = {-irT{q+l)T{-q)- 


\/ 2 'kx 


(l + 0(x-’)) . 


^It has the asymptotic behaviour: 

'd(t|a) ~ for t oo ] d(<|a) ~ for t ^ 0 . 


(B.16) 


(B.17) 


(B.18) 


(B.14) 
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where the last expression has a meromorphic extension with the same general properties 
as those of (^{s,d). We remark that the reflection formula does not hold in the general 
case a 7 ^ 0. 

C Properties of dimensional continuation of sums 

We present here some basic properties of the dimensionally regularized series as defined 
by (|l^. They are the analogues of those discussed in |p, valid for dimensionally regular¬ 
ized integrals. 

Property 1. 


neZP mez? neZP+9 


Proof. Consider a Mellin-transformable function f{x), assuming for simplicity that 
/(x) —> 0 for X —> 0 and /(x) ~ x~^, p > 0 for x —> oo, such that its series can be 
defined without IR subtractions, which are irrelevant to this discussion. Let us denote by 
M[f, s; y] the Mellin transform of /(x -|- y) with respect to the variable x: 

poo 

M[f,s-y]= dx/(x-Fy)x*"^ . 

Jo 

Transforming also on y and performing a change of variables we get 


poo poo pi 

M[f,s;w] = / dyM[f,s-,y]y^-^= / dtf{t)P+^-^ dvv^-\l-vr-^ 
Jo Jo Jo 


r(s)rM 
r(s -i- w) 


(C.2) 


J^[f,s + w] , 


with Re s, Re tc > 0, Re s -|- Re tc < p for these expressions to make sense. Applying the 
definition of regularized sums and using (|C.2|), it follows that 


m£ZP neZi 


1 

(27r 



ds C{w,p)C{s,q) 


ns)r{w) 
r(s -i- w) 


J^[f,s + w] 


By making a change of variable in the double complex integral and taking into account 
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eqs. (p.8D and (B.£), we get 


' Yl' = ( 2 ni )2 [ [ dwT{w)C{w,p)T{u-w)C{u-w,q) 


m£ljP nGlj^ 


^ / duM[f,u] [-C{u,p)-C{u,q) + C{u,p + q)] 
27ri Jr^ 

- Y' ~ Y' Y ’ 

m&P n&Zi neXP+i 


(C.3) 


where r„ is a contour laying in the half-plane Re w < p. This is exactly the equality we 
are looking for. 


Property 2. 


J d^p Y' ^Y'j ^^p f^p^' 


n 


(C.4) 


Proof. It is convenient to introduce an auxiliary function 

/(a,p2,n2) = /(p2,n2)e-“(?’^+-^) , 


for which the right- and left-hand sides of ( |C.4| ) become 

l.h.s. : ^ J (|C^ ) 

^ fia,x,y) . 


r.h.s. 


(|C.5|R) 


Thanks to the exponential factor, the series and the integral in both expressions have been 
defined without recourse to subtractions, by simply choosing a large enough value for (the 
real part of) D and d to have IR convergence. Also, the contour T = {Re s = c; c > d/2} 
can be defined to be the same in (C.5L) and (|C.5[ R) by fixing c sufficiently large. Then, 
property 2 for the auxiliary function follows trivially from exchanging the integrals over 
X and y and ( C.5 L), ( C.5| R) are both equal to the same function I{D, d, a) analytic in its 
variable. The analytical continuation of I{D,d,a) down to smaller D and d is still given 
by d^), Ol), but now with subtractions made. If we take (the real part of) D and 
d small enough to have UV convergence even without the exponential factor, we can put 
a = 0 and (p.4|) follows. 


Property 3. 


Y'Y'= Y'Y' 

neZP m&zi meZi n£ZP 
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Proof. We give the proof only for the simpler case in which the function 
does not depend on the product {m ■ n). Proceeding as in the previous case, one can 
introduce the auxiliary function 

/(a, , 


for which the right- and left-hand side of (p.6|) become 



i ds Cis,p) 1 

Vs 

r 

^ dx / 

0 J 

poo 

du C{u,q) 1 
J 

r 

poo 

' dy f{a,x,y) (|CU|L) 

0 

poo 

1 • 11 . O . . , X J 

{2mY J 

r.h.s. : ^ j 

(27rf)2 J 

1 du C{u,q) 
r„ J 


1 ds Cis,p) j 

Ps tJ 

I dx x^~^ f{a,x,y) (|C.7|R) 
0 


Again, both series have been defined without recourse to subtraction, by choosing a large 
enough value for Re s and Re u along the contours P* = {Re s = c; c > p/2}, P. = 
(Re u = c]c > q/2}. Property 3 follows trivially for the auxiliary function by simply 
exchanging the various integrals and we recover o by continuing both sides (E3L), 
( C.7 R) analytically to o = 0. 
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